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Abstract 



We investigate the effects produced by tfie tliree-momentum scale Y^mAqcD in the 
strong coupling regime of heavy quarkonium. We compute the leading non- vanishing 
contributions due to this scale to the masses and inclusive decay widths. We find 
that they may provide leading corrections to the S-wave decay widths but only 
subleading corrections to the masses. 
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1 Introduction 



Heavy quarkonium is characterized by the small relative velocity v of the heavy quarks in 
their centre-of-mass frame. This small parameter produces a hierarchy of widely separated 
scales once multiplied by the mass m of the heavy particle: m (hard), mv (soft), mv'^ (ultra- 
soft), .... In general, we have E ~ mv^ <S p ~ mv <S m, where E is the binding energy and 
p the relative three momentum. 

It is usually believed that for most of the heavy quarkonium states a weak coupling analysis 
is not reliable. However, one can still exploit the hierarchy of scales in the problem [1]. It 
was argued in [2,3] that in the particular case Aqcd ^ mf^, which we will be concerned 
with in this letter, it is possible to encode all the relevant information of QCD in an effective 
Schrodinger-like description of these systems. The problem then reduces to calculating the 
potentials from QCD. It has been shown in [3] how to systematically calculate the potentials 
within a 1/m expansion (see [4] for earlier calculations). 

Once the methodology to compute the potentials within a 1/m expansion has been developed, 
the next question appears naturally: at which extent one can compute the full potential 
within a 1/m expansion in the case Aqcd ^ mv"^.^ We tackle this issue in this paper. 
We will see that, indeed, new non-analytical terms arise due to the three momentum scale 
Y^mAqcD- These terms can be incorporated into local potentials (5^(r) and derivatives of it) 
and scale as half- integer powers of 1/m. Moreover, we show that it is possible to factorize 
these effects in a model independent way and compute them within a systematic expansion 
in some small parameters. 

As mentioned before, these terms are due to the existence of degrees of freedom, namely the 
quark-antiquark pair, with relative three momentum of order y^mAqcD- The on-shell energy 
of these degrees of freedom is of O(Aqcd), i-e. the same energy scale that is integrated out 
when computing the standard 1/m potentials, which corresponds to integrating out (off- 
shell) quark-antiquark pairs of three momcntTim of order Aqcd- Therefore, in principle, both 
degrees of freedom should be integrated out at the same time. 

In this letter, under the general condition Aqcd ^ mv"^, we will perform the analysis in two 
possible cases: 

1) in Sec. 2 we will consider the particular case mv ^ Aqcd; 

2) in Sec. 3 the general case Aqcd < "mv. 



Note that the scale y mAqcD fulfils ymAQCD ^ n^v and y'mAqcD S> Aqcd- From the last 

^ In fact, there is at least one example where powers of ^/m arise upon integrating out some 
non-relativistic degrees of freedom [5] . 
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inequality it follows that at this scale we always are in the weak coupling regime. 



2 Case mv Aqcd 



In the case mv ^ Aqcd, all quarks and gluons with energy much larger than Aqcd (in 



particular gluons with energy and momentum of order ^ttiAqcd and mv) may be integrated 
out from NRQCD using weak coupling techniques. This leads to the EFT called pNRQCD' in 
[6,7] (formerly called pNRQCD in [8,2]). This EFT contains, as explicit degrees of freedom, 
gluons with energy and momentum smaller than mv and quarks with energy smaller than 
mv and momentum smaller than m. Quarks may be arranged in quark-antiquark singlet 
S = Slc/VK and octet O = l/^O^'T" fields {Tp = 1/2). The Lagrangian of pNRQCD' 
then reads (R is the centre-of-mass coordinate and r the relative coordinate) [2]: 

I^pNRQCD' = Jd'R J d^'r (^Tr jst (^^o - /i.) S + {tDo -ho)0^ 

+Tr{0^r-^ES + S^r-EO} + iTr{0'^r-c/EO + O^Or-^E}j + J d^RLg, (1) 

where hg — — V^/m + Vs, ho — — V^/m + K, and jCg stays for the Lagrangian density 
of gluons and light quarks. The potentials V = {Vs,Vo} contain real and imaginary parts. 
The real part, which at leading order is the Coulomb potential V'^'^\ has been calculated by 
different authors over the past years [9]. The imaginary part has been calculated in [6,7]. 
It consists of local potentials {6'^{r) and derivatives of it). The imaginary coefficients come 
from the imaginary parts of the four-fermion matching coefficients of NRQCD [1] . 

The next energy scale to be integrated out is Aqcd- This means integrating out all quarks and 
gluons of energy or kinetic energy of order Aqcd- The contributions due to (off shell) heavy 
quarks of energy ~ Aqcd and three momentum of order mv or smaller (i.e. of order Aqcd) 
are easily singled out by performing an expansion of the incoming and outgoing bound-state 
energies hs and ho over Aqcd in the matching calculation. This ensures that the quark kinetic 
energy is much smaller than Aqcd and, therefore, that the quark three-momenta are much 



smaller than JmAqcB- This expansion only produces terms that are analytical in 1/m [6,7]. 



The contributions due to heavy quarks of three momentum of order y'mAQCD may be ob- 
tained as follows. We spht the singlet and octet fields of the pNRQCD' Lagrangian into two 
fields: 

s = Sp + s,,, o« = o; + o:„ (2) 

where the semi-hard fields Ssh and O^j^ are associated to three-momentum fluctuations of 
O ^^tbAqcd) and the potential fields Sp and to three-momentum fiuctuations of 0{mv). 
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The potentials arc labeled according to the relative momenta that they connect: V = V^'^ + 
y p,sh y sh,p ^ ysh,sh_ rpj^g typical thrcc-momcntum transfer in V^^^^, V^'P and V^'^'^'^ is 

^mAqcD (> mv). 

The pNRQCD' Lagrangian then reads 

LpNRQCD' — Lg + L^pNRQCD' + ^pNRQCD' + ^mixing- (3) 



The expressions for L^J^j^rqcd' ^pNrqcd' identical to the pNRQCD' Lagrangian 

except for the changes S, 0", K ^.ft, O^/., Vf'''', V„'^''^ and 5, O", K, K ^ S'p, 
Op, Vf'P, V^'P respectively. Recall that the gluons left dynamical are of O(Aqcd) and that 
analytical terms in r do not mix semi-hard and potential fields. Therefore, the multipole 
expansion in (1) is an expansion with respect to either the scale r ~ 1/ ^JmAqcD in -^pAfRQCD' 
or the scale r ~ 1/mv in L^p^rqcd'- 

Throughout the paper we will also assume that 



mAqcD > masUmAqcD) , (4) 



which implies that the Coulomb potentials in V^''"^^, V^^'"'^ and V^^'''^^ can be expanded about 
the kinetic energy and no Coulomb resummation is needed. This is not so for V^'^. 

The leading contribution to the real part of Lmixing comes from the mixing of Ssh with Sp 
and O^^ with due to the Coulomb potential. As an example, consider the real part of the 
singlet-mixing term due to the static Coulomb potential. It is given by 



-^mixing 



Singlet 



- j dPR J d\ Sl{K,r)VP{r) Ssh{'R,r) + R.C. (5) 

--Jd'Rj d'p J dV^J(R, P) K^'Kp - P') S,,{R, p') + H.c. 

^-Jd^R J d^r (^Sl{R, 0) + r • Vr5't(R, r) + • • • j V}""^ (r) S^hCR, r) + H.c. . 

In order to avoid a cumbersome notation we have dropped the upper-index p. sh from Vj'^\r). 
In fact, any potential between fields labeled by a, 6 = p,sh always has upper- indices a,b. 
Hence, dropping the upper-indices shall not lead to ambiguities. In the second line of Eq. 
(5), a Fourier transform of all the fields has been performed, and in the third one, we have 
expanded around p ~ in the potential, since, by definition, p ~ mv ^ p' ~ y^mAqcD- 
Doing so in the loops that will appear in the matching computation guarantees that only the 



scale ymAqcD is integrated out. Alternatively, one may consider S'j!j(R, r) slowly varying in 
r and multipole expand it about r = 0, which brings us directly from the first to the last 
line of Eq. (5). At the order of interest we have Vj^^ = —Cj a^jr and = ttg (^tbAqcd)- 
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Analogous results hold for the real part of the octet-mixing term due to the static Coulomb 
potential: 



R-G Ljjiixing (6) 
Octet 

^-Jd'R J c^VTr|(|ot(R,0) + r- VrOj(R,r) + • • • j V;W(r) O^^lR, r) + H.cj , 
where the trace is over the colour indices, the mixing potential is V^*^''-* = l/{2Nc)as/r and 

The leading contribution to the imaginary part of L^ixing can be immediately read off from 
the imaginary delta- type potentials calculated in [7]: 



ImL„,,,,, = - I d?R I d'rTr |sI,(R, 0) ^<5=^(r) Sp(R, 0) + H.c.j 

- J d'R J d\Tr jo^jR, 0) ^5^v) 0,(R, 0) + H.c.j , (7) 



where 



+4Im/EMC-5o) -2 8^ (Ihi/emC-^o) -Im/EM('-5i)) ), (8) 



Ko^-'^I^Um fsCSo) - 2 (im fsCSo) - Im fsi'S,)) J . (9) 

The matching coefficients / are the matching coefficients of the four-fermion operators in 
NRQCD and may be read off from Ref . [1] . 



2.1 Matching 



The next step is to integrate out from pNRQCD' all fluctuations that appear at the energy 
scale Aqcd- These are light quarks and gluons of energy or three momentum of order Aqcd, 



and singlet and octet fields of energy of order Aqcd or three momentum of order y'mAQCD- 
We will be left with pNRQCD, where only a singlet field describing a quark-antiquark pair 
of energy mv^ and relative three- momentum mv is dynamical ^ : 



We ignore pseudo-Goldstone bosons (pions), which, in principle, should also be included. 
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LpNRQCD' LpNRQCD — LI^NRqcd + ^l^NRQCD ) (10) 

L%QCD = Jd'Rj d'r (^td^ - ^ - V^' - 5^^/™) S, (11) 

L'Jn^cd = -jd'Rj d'r SV'^^ S. (12) 

L^Jnrqcd is defined as the part of the pNRQCD Lagrangian obtained by integrating out 
quarks and gluons of energy and three-momentum of order Aqcd in ^pNRQCD' only. It is 
analytical in 1/m and has been considered before in [2,7]. Here we will calculate the leading 
part of L^Jnrqcdi which is defined as the part of the pNRQCD Lagrangian obtained by 



integrating out quark-antiquark pairs of three-momentum y mAqco in Lpf^RQco' in addition 
to the above degrees of freedom. In general, it is non analytical in 1/m, and, at leading order, 
it consists of a new local (delta-type) potential. 






Fig. 1. The four diagrams of pNRQCiy contributing to SV^I"^ at leading (non-vanishing) order in 
the multipole expansion. Open and full circles indicate octet and singlet potential insertions coming 
from the mixing terms respectively. These are treated according to Eq. (5). The upper-scripts P 
and SH on a propagator indicate that the propagating fields are of the potential and semi-hard type 
respectively. The circle with a cross indicates the vertex S^^ r • E O (or Hermitian conjugate ), where 
the quark fields are both either potential or semi-hard. The gluon line stands for non-perturbative 
multi-gluon exchanges. 



The matching condition for the full 5V 
in the multipole expansion is 



1/m _|_ ^yi/ -J^ leading (non- vanishing) order 



-5V- 



E - K K. 



— / dt (vac|ir • gE{t) e 



-i{ho-E)t 



ir ■ 5(E(0)|vac) 



E-K 



(13) 



The above matching equation should be understood (even if written at the operator level) 
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with incoming (outcoming) momenta and energy E of 0{mv) and 0{mv^) respectively. The 

typical size of the time variable in the integral is given by the vacuum expectation value of 
the chromoelectric correlator and hence t ~ 1/Aqcd- The separation between potential and 
semi-hard relative three-momenta discussed above can be easily implemented in the rhs of 
Eq. (13) by expanding the Hamiltonians hg^o in K^o'^ V^'^^. The zeroth order term in 
this expansion gives 5V^/'^ and has been calculated in [2,7]. The V'^''^ potential cannot be 
expanded in the potential region. The size of the three-momenta in the semi-hard regions 



is of 0{^mA.qp,T>)- Several approximations apply: (i) E — ho,s ~ — p^/m ~ Aqcd in the 

semi-hard regions, (ii) {ho — E)t ~ + V^'^ — E)t ~ mf^/AqcD ^ 1 in the potential 
regions and (iii) we can expand the potential three-momenta with respect to the semi-hard 
ones in V^^'^ and V^'^^. 

The leading contributions to 5V^/^ have been depicted in Fig. 1. Fig. la corresponds to 

oo t 

^yVV^ = 1_ j dt Jdt' (vac|ir • gE{t) (-iVP''^') e~'^''ir ■ ^E(0)|vac)^l//'^'^', (14) 

'^0 m 

Fig. lb corresponds to 

. oo t 

Sv^^ = l/p.*/^—-! Jdt J dt' (vac|ir • gE{t) e"^- {-iVf'^) ir ■ yE(0)|vac), (15) 



m 



Fig. Ic corresponds to 



oo 



Syi/V^=j_ j dt j dt' J dt" (vac|ir • gE{t) 

■^0 

X (-iVP''^) e"* (-iV,''''P) ir ■ ^E(0)|vac), (16) 
and, finally. Fig. Id corresponds to 

. oo 

5yi/v^ = f dt (vac|ir • gE{t) e-'^Ur ■ gE{0)\vcic)—Vf^P. (17) 



m 



The potential 6V contains a real and an imaginary part. The real part contributes to the 
heavy quarkonium spectrum, the imaginary one to the inclusive decay width. 
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2.2 Corrections to the spectrum 



The four diagrams that give the leading contribution to Re 5V^/ ^ are obtained from those of 
Fig. 1 by substituting Vg^o K^,o^ where V^^] are the Coulomb singlet and octet potentials. 
They give: 

= (2C; + CA?^^^7:al£f/, ^ , (18) 

where in the first equality we have used the definition of that one may find in Ref. [7] and 
in the last equality we have written the chromoelectric correlator in Euclidean space (traces 
as well as suitable Schwinger lines connecting the gluon fields are understood): 

1 °° 

^" " TV / ^^^"(vac|5E(r) ■5E(0)|vac)e. (19) 







Eq. (18) gives a contribution to the energy of O I mv as 



mA 



QCD^ 



2.3 Corrections to the decay width 



The four diagrams that give the leading contribution to ImSV^^"^ are shown in Fig. 2. 
These can be derived from the diagrams of Fig. 1 by replacing one of the potentials by a 
Coulomb potential and the second potential with the imaginary delta potential of Eq. (7). 
The graph with two potentials inside the gluonic loop as well as graphs involving the octet 
delta potential (~ KoS^{r)/m'^) do not contribute to ImSV as a delta potential (although 
they do as derivatives of a delta potential, which are subleading). We obtain 

Im SV^'^ = -z^/^^ {2Cf + Ca)^ a, £5/2 ^ 
45 -v/TT m^/^ 



where in the last equality we have written the chromoelectric correlator in Euclidean space. 

The above correction gives a contribution to the decay width suppressed by a factor of 
O I ^'^^ mas _ I ^j^]^ respect to the leading contribution. 
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V K K V 

OS so 





Fig. 2. The four diagrams of pNEQCL/ contributing to Im(5F^/"^ at the leading (non-vanishing) 
order in the multipole expansion. The full box indicates the insertion of a delta-type potential pro- 
portional to Kg. All other symbols are as in Fig. 1 with Vs,o ^ Vj^j . 

A similar analysis can be done for the P-wave decays. The leading effect would be in that case 



at least 0{mas/ ymAqcD) suppressed with respect to the leading contribution computed in 
[6]. 



3 Case Aqcd < mv 



Here we will follow the same procedure as in the previous section. In this case, however, 
the starting point is the NRQCD Lagrangian. Wc split the quark (antiquark) field into two: 



a semi-hard field for the (three-momentum) fluctuations of O(ymAQCD), '^sh (Xsh): and a 
potential field for the (three-momentum) fluctuations of 0{mv), ipp {Xp}'- 

= V'p + i^sh, X = Xp + Xsh- (21) 



The NRQCD Lagrangian then reads 

sh j_ TP 

'NRQCD "T ^ NRQCD 



L NRQCD — L^NRQCD + ^ NRQCD + ^mixing + Lg. (22) 



The Lagrangians Lf^p^QQj-, and L^j^irqcd identical to the NRQCD Lagrangian expressed in 
terms of semi-hard and potential flelds respectively. The quantity Lg is the QCD Lagrangian 
for gluons and light quarks. For L^j^hqcd can use weak coupling techniques. Therefore, 
we can construct a pNRQCD' Lagrangian for it, once gluons and quarks of energy or three 
momentum of 0{^m Aqcd) have been integrated out and transformed into potentials: 

^ NRQCD ~^ ^pNRQCD'- (23) 
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If we further project to the quark-antiquark sector, the Lagrangian L^pf^RQCW ^'^^ formally 
read equal to Eq. (1). The multipole expanded gluons in L^'j^j^QQjj, have (four) momentum 

much smaller than ^JmK^cD■ We note that we cannot do the same for LFj^rqcd since at 
scales of O(Aqcd) we can neither use weak coupling techniques nor the multipole expansion. 




Fig. 3. The Coulomb-exchange graph contributing to the leading mixing interaction between 
semi-hard and potential fields. 

We consider now Lmixing- We will assume, as in Sec. 2, that the condition (4) holds. This 
will allow us to treat the Coulomb potential as a perturbation at the semi-hard scale. The 
leading order contribution to the real part of Lmixing comes from the one-Coulomb exchange 
graph (see Fig. 3): 



ReL^Xrng = - j R J T Tt {j\R) vj'>\v) S.,(R, r) } + H.c. 

-J d?R j ciVTr{jt(R)i;(0)(r)O,;,(R,r)} + H.c. , (24) 
jt(R)^Xp(R)V'J(R)- (25) 

The potentials V^^^ and V^;^^ are perturbative: V^'^'> = —Cfag/r and Vj°^ = l/{2Nc) ag/r. 



The coupling constant is calculated at the semi-hard scale wmAqco- Besides the above term 



wc need to consider also the next-to-leading term in the mv/ ^JthAqcd expansion. It is given 

by 

ReLS.,„, = - / d'R J ciVTV{jt(R).ryio)(r)S,,(R,r)} + H.c. 

-Jd'R ydVTV{jt(R).rV;(°)(r)0,,(R,r)} + H.c. , (26) 

jt(R)^Xp(R)^^^'J(R)- (27) 

A practical way to obtain ReL^l^j^^^ is by expanding the Coulomb potential in Fig. 3 at 
higher order in p/p' and promoting the conventional derivatives acting on the potential 

fields to covariant ones. A proper tree-level matching in coordinate space can be done using 
the field redefinitions of Ref. [8] for the semi-hard fields projected to the two particle sector 
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and multipole expanding the potential fields. The leading contribution to the imaginary part 
of L^ixing is analogous to the one given by Eq. (7): 



ImL(^Lnp = - / d'R J d\TTl^Sl{K,0)^S\r) J(R)| + H.c. 

-Jd'Rj d'r Tr[oUR, 0) ^5^(r) J(R)} + H.c. . (28) 



Note that the potential fields always appear as local currents in Lmixing- Finally, the effective 
field theory that we obtain, at the order of interest, is given by 

LnRQCD' — Lp%iiQ(jD, + L^RQCD + ^Uixing + ^mixing + I™ ^mixing + -^5 ) (29) 

where Lg contains now gluons and light quarks of energy and momentum much smaller than 
^Jm Aqcd- 

3.1 Matching 



As in section (2.1), we now want to integrate out degrees of freedom of O(Aqcd)- We will 
be left with an EFT, pNRQCD, where only a singlet field describing a quark- antiquark pair 
of energy mv'^ and relative three momentum mv is dynamical: 

Lnrqcd' LpNRQCD = L^Jnrqcd + L^Jnrqcd- (30) 

The quantity L^Jnrqcd obtained by integrating out quarks and gluons of energy and 
three momentum of order Aqcd in L^nrqcd- i^ analytical in 1/m and has been considered 
before in [3,6,7]. Here we will calculate the leading part of L^n^cd^ which, in general, is 
non analytical in 1/m. It involves the integration from NRQCD of quark- antiquark pairs 
of three momentum y^mAqcD- The Lagrangian L^Jnrqcd consist, at leading order, of a 
new local (delta-type) potential that we name 5V^/"^: 



LpNRQCD 



^- j d^nj dh S"^ 5V^'^ S . (31) 



The matching calculation for 5V^/^ is analogous to the computation of the previous section 
supplemented with the technology developed in Refs. [3,6,7]. The leading contribution is given 
by the four diagrams shown in Fig. 4. Fig. 4a corresponds to (according to the notation of 
Ref. [3]) 

-X2)53(xi -x;)<53(x2 -x'J = j d'R (°)(0;xi,X2| J'(R) 
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Fig. 4. The four diagrams of NRQCL/ contributing to SV^I"^ at leading order. Full circles indicate 
singlet potential insertions coming from L^mixing- -^^^ upper-scripts P and SH on a propagator 
indicate that the propagating fields are of the potential and semi-hard type respectively. The circle 
with a cross indicates the vertex S^^r • EOgh (or Hermitian conjugate). The box with a cross 

indicates the vertex Ke L^^-^-^g. The gluon line symbolizes multigluon non-perturbative exchanges. 



(p = 0|r'i;(°)(r 



Fig. 4b corresponds to 
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— El — }J _p1 



V;(o)( 



r) |p = 0) 



<7E^(R) J(R)|0;x;,x'2)(°), (32) 



5yi/v^(xi - X2)53(xi - x;)53^X2 - x'2) = j £R (°)(0; xi, X2I yE^R) J(R) 

J^(R) |Q;x;,x^)(°), 



(p = 0| Vi°)(r) -^-j^ r= Vf\r) |p = 0) 

m m 

Fig. 4c corresponds to 

(5yV^(xi - X2)53(xi - xi)(5^(x2 - x'2) = J d?R (°)(0;xi,X2| J^(R) 

J^(R) |0;x;,x'2)(°). 



(33) 



(p = 0| V^'\t) J v= V^'\t) Ip = 0) 
Finally, Fig. 4d gives 



m 



(34) 



W^'^i^i - X2)5^(xi - x;)5=^(x2 - x'2) = j d^R (°)(0;xi,X2| t/E^(R) J(R) 



X 



(p = o|v;(°)(r) 



m m m 
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x^E^(R) J(R)|0;x;,x^)(°). (35) 

In the above equations H stands for the NRQCD Hamiltonian in the static hmit and |0; r)'^'^) 
is the gluonic piece of the ground state of NRQCD in the static hmit. We refer to [3,7] for 
further details. As in Sec. 2, several approximations apply: (i) /io,s ~ —p^/m H Aqcd 
in the semi-hard regions, (ii) whereas hg^s <S -f^ ~ Aqcd in the potential regions. Moreover, 
we can expand the incoming (outcoming) three-momenta with respect to the semi-hard ones 
in Vfll'P and VP'^'^. 

By summing up all the contributions, we obtain the same result as in Sec. 2: 

Re SV'/^ = {2Cf + C^f^^^^ ^ . (36) 



This is not a coincidence. Note first that the diagram in Fig. 4d is identical to the one in 
Fig. Id. The remaining diagrams in Fig. 4 also have a mapping to the corresponding ones 
of Fig. 1, if we substitute the square box in the former by a round box linked to an open 

circle through an octet propagator. This mapping can be made rigorous from the following 
equality (where {l^^)*-"-*} is the gluonic term of a complete set of eigenstates of the static 
NRQCD Hamiltonian, and E^^ the corresponding eigenvalues [3,7]): 



I d^R (°H0;xi,X2| J(R) = (°)(0| D,,, 53(xi - X2) 
= ^(°)(0|D.,53(xi-X2)|n)W + 

= 53(xi - x^)-^ i:(vac| 5E(xO-75y^-|n)(°) W(n| + . . . , (37) 

where in the last line we have also made use of the fact that in the limit Xi — X2 — ^ we 
have |0)*^°-' lc|vac)/-\/]Vc. The neglected terms, gencrically denoted with dots, do not give 
delta-type contributions to the potentials. From Eq. (37) it follows that the calculation of 
the diagrams of Fig. 4a, Fig. 4b and Fig. 4c reduces to that one of the diagrams of Fig. la. 
Fig. lb and Fig. Ic respectively. Similarly, for the imaginary part of bV^I"^ the relevant 
diagrams reduce to those calculated in Sec. 2.3 and shown in Fig. 2. It reads: 

Im W'l^ = (2C, + Ca) a. ^ . (38) 



Here, as well, an analysis for the P-wave decays could be done. We can easily estimate 
that the leading effects would be at least O {mas/ \J m^qcn) suppressed with respect to the 
contributions computed in [6]. 
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4 Conclusions 



For heavy quarkonium systems in the strong-coupUng regime (Aqcd ^ mv'^), the corrections 
to the static QCD potential in the Schrodinger equation have so far been calculated within 
a 1/m expansion. We have shown here in a quantitative manner that they are not the only 
contributions to the full potential and have computed the leading non-analytical corrections 
in 1/m. 

Our findings can be summarized in the following corrections to the energy levels and the 
(S'-wave matrix elements and decay widths (the symbols V and P stand for the vector and 
pseudoscalar S'-wave heavy quarkonium respectively, n is the principal quantum number): 



SE = {2Cf + CaT 



3r(9/2) 



2 cE |-^nz(0)|% 
«s ^7/2 ^t:^^, 



3/2 



{VQ{nS)\0,{'S,)\VQ{nS))^CA- 



2n 
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^ 3r(7/2) mV2 ^^\m 



(39) 
, (40) 



{PQ{nS)\0,{'S^)\PQ{nS))^CA 



I T^P 



27r 



\ 3r(7/2) mV2 ^"^Vm 



, (41) 



r{VQ{nS) ^ LH) = 



r{PQ{nS)^LH)- 



TT 

X 



CA|i?ro(0)P 



TT 

X 



+ 



, (42) 



, (43) 



T{VQ{nS)^e^e-)- 



TT 



X 



Im/ee('-Si) 1 + 



4(2Cj + Ca) «sg5l2 ^ ^/l 
777,1/2 



3r(7/2) 



r(PQ(n5)^77) = 



^A|i?^o(0)P 



TT 
X 



+ 



,(44) 



,(45) 
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where 0{l/m) stands for corrections (which may be of the same size) that can be computed 
within the 1/m expansion (see [7]) and for higher-order corrections. 

Let us comment on the size of the new corrections. For the spectrum they are always smaller 
than mv^ and therefore subleading with respect those calculated in [3] . For the S-wme decay 
widths their relative size with respect the corrections computed in [7] depends on the size of 

asi^mAqc]:,). Under some circumstances, for instance ~ v, the contributions calculated 
here are the dominant ones. In any case, the above results fulfil the same factorization 
properties as those obtained in [7]. As a consequence, equations like those given in Sec. Vll 
of Ref. [7] still hold. Let us also note that the same non-perturbative correlator appears in 
both electromagnetic and hadronic decays. 

In this paper we have assumed that the scale y^mAgcD is much larger than mag. Otherwise 
we are not allowed to treat the Coulomb potential as a perturbation at that scale. This may 
not be the case for the T system where one seems to be in the situation Aqcd ~ mctg, 
which implies ^ mAgcD ~ mccg. In this case, one should integrate out the three- momentum 

scale mas at the same time as the scale ^JmAqcD- The calculations presented here should 
be modified by using the full Coulomb propagators instead of the free ones in the semi-hard 
regions. In addition extra contributions may arise, which are only due to the three-momentum 
scale mas- We do not deal with this issue in this paper, which, however, deserves further 
studies. 
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